Region of fidelities for a 1 — >■ N universal qubit quantum doner 
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We analyze a region of fidelities for qubit which is obtained after an application of a 1 — > N 
universal quantum doner. We express the allowed region for fidelities in terms of overlaps of pure 
states with irreps of S n in = N + 1) showing that the pure states can be taken with real coefficients 
only. Subsequently the case n = 4, corresponding to a 1 — > 3 doner is studied in more details as an 
illustrative example. To obtain the main result, we make a convex hull of possible ranges of fidelities 
related to a given irrep. The formalism allows to construct the state giving rise to a given N-tuple of 
fidelities. 
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I. INTRODUCTION 

A. Historical overview: the no-cloning theorem 

In 1982 Nick Herbert published a paper about faster 
than light communication, based on quantum correla- 
tions. He called his project FLASH, an acronym for 
'First Laser- Amplified Superluminal Hookup' |JJ. Of 
course, the idea was incorrect but it was the source for 
Wootters and Zurek [2 J and independently Dieks J3j to 
establish the so-called no-cloning theorem. The proof 
(we follow the version of the proof as in (4|, |5J, ||6l ) 
is the following; consider a cloning unitary operation 
U c i that is able to clone every given qubit. Before the 
copying process, we have states: 

\A) = \x)\0)\M), 

\B) = |x')|0)|M), (1) 

where \x), \x) denote 'text' to copy for each state, |0) 
represents a 'blank card' and \M) is a cloning machine 
state (M denotes the cloning machine). Acting U c / on 
([I}, we should get: 

u cl \A) = I*)I*)|m'), 

U cl \B) = \x')\x')\M"), (2) 

where |M ), |M } represent new states of the machine. 
Now let us consider a scalar product in both cases: 

(A\B) = (x\x), (3) 

since (0|0) = (M\M) = 1. On the other hand, for the 
second case, we have: 

(U cl A\U cl B) = ((x\x')) 2 (m'\m"). (4) 



email: studzinski .m. g@gmail . corn! tel: (+48 58) 551-20-34 



Comparing |3) and Q, we get that in the case of 
< \{x\x )| < 1 scalar products do not match each 
other (but from the property of the unitary action, they 
should be the same). So, we have contradiction and it 
proves that one is not able to clone initially unknown 
quantum states. In general, the no-cloning theorem 
comes from the linearity of quantum mechanics |2|, |3|. 
Note that the no-cloning theorem could be generalized 
||7| or even strengthened [8j. 



B. Beyond the no-cloning theorem: origins of quantum 
cloning 

Although, the no-cloning theorem is fundamental for 
quantum physics, it is not of great use in practice. The 
no-cloning theorem states that one is not able to copy an 
arbitrary quantum state. However, it is well known that 
performing ideal operations in physics and especially 
in quantum physics is impossible. On the other hand, 
it is obvious that one can copy, though perhaps with a 
very bad quality. Thus it is crucial to know the ultimate 
bounds for the quality of copying. 

Several years after Wootters and Zurek paper has 
been published, namely in 1996, Hillery and Buzek 
published a paper called 'Quantum copying: beyond 
the no-cloning theorem' J9J- It was the first time, 
when the above question regarding imperfect cloning 
has been formulated. Subsequently, the subject was a 
matter of wide research (see |7| and references therein 
for a comprehensive review). Here, we would like to 
point only a few (regarding the cloning of finite quan- 
tum states only). Quite soon, the Buzek-Hillery 1 — s> 2 
(qubits) Quantum Cloning Machine (QCM) (for all for- 
mal definitions of quantum cloning machines, we refer 
to [7J), was generalized to the case N\ — > N2', first for 
qubits by Brufi et al. in (lOl , and by Gisin and Mas- 
sar in [llj, and then for arbitrary-dimensional states 
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by Werner in [12], and Keyl and Werner in |13|. Thus 
the family of symmetric Universal Quantum Cloning 
Machines (UQCM) is, at present, well known. What 
is more, the asymmetric UQCM were also very heav- 
ily studied, let us mention here works of Braunstein 
et al. Ql, Cerf 03, Fiurasek et al. QSI, and Iblis- 
dir et al. [17], |18|. To emphasize the relevance of all 
these works, let us mention that formalism of asym- 
metric cloning machines is important in the context of 
quantum cryptography - it can be used in studies of 
relations between the eavesdropper's information gain 
and the noise in the channel. Later, efforts were made 
to unify these two kinds of QCM [19]. Of course, a lot 
of questions and problems still need answers, for in- 
stance, optimal state-dependent QCM (as an example, 
see 11201 or ET] - unfortunately, not many results are 
known for this kind of QCM) or optimal asymmetric 
QCM. There is also another 'gap' in quantum cloning: 
interestingly, up to our best knowledge, one is lacking 
a general result on an admissible region of fidelities (in 
general not the optimal one) for universal asymmetric 
1 — » N quantum cloning machines; in our work, we 
want to make progress in this direction, providing an 
answer to this particular problem in the case of qubits. 
Let us mention briefly that it has been studied partially, 
for example in [ 16 1 (and partially in [17]), where expres- 
sions for an optimal fidelities region for an asymmetric 
1 — >■ 3 quantum doner were obtained (and some partial 
results for a 1 — > N doner). 

In this Letter, we shall consider a 1 — > N universal 
quantum cloning machine (quantum doner) for qubits. 
We first point out that this problem could be related 
to singlet states by recalling a relation between fideli- 
ties and singlet states. Using this fact, it turns out that 
by the application of Schur - Weyl duality, our prob- 
lem could be quite easily solved and it leads to plots of 
ranges of fidelities for different irreducible representa- 
tions of the symmetric group S n (where n = N + 1). 
The method of irreps (irreducible representations) is 
'nice' here, because it makes calculations a lot easier. 
It allows us to decompose our initial Hilbert space into 
blocks of smaller dimensions, connected to a given par- 
tition A (A's label irreps of a symmetric group), and 
moreover, in our calculations, we can restrict ourself to 
a pure state only, linked to a given block. After taking 
the convex hull of the figures corresponding to all par- 
titions, we can obtain the possible range of fidelities. In 
our case-study example of the 1 — > 3 doner, we check 
that our results are consistent with the existing meth- 
ods for 1) symmetric doners (for example, Keyl and 
Werner work [13]), namely we obtain that in the case of 
optimal, symmetric UQCM, result of all three fidelities 
equals | is obtained, 2) an optimal asymmetric 1 — > 3 
doner obtained in 1 16 1 . Finally, as a direct application of 
our method, we study one particular example, namely 
the case when one has max (F\ + F3 = 2F2). 

This work is organized as follows. In Section [IT] we 



start with showing that the quantum cloning problem 
is equivalent to the entanglement sharing picture, so 
that the cloning machine is equivalently represented by 
a multipartite quantum state. Then, we formulate prob- 
lem that we want to address in this Letter, namely, to 
determine the region of allowed fidelities for our QCM. 
What is more (and crucial for us), we show a transfor- 
mation between Bell's states that allows to use Schur - 
Weyl duality. Section [III] is devoted to the group rep- 
resentation theory, in particular, Schur - Weyl decom- 
position. We present the basic formalism that leads to 
Schur - Weyl duality. At the beginning, we introduce 
SWAP operators and then we show how they lead to 
Schur - Weyl decomposition. At the end of this sec- 
tion, we show how Schur - Weyl duality is connected 
with Young diagrams formalism. This section is strictly 
mathematically oriented, but as it will be shown in the 
next section, this formalism allows us to simplify our 
problem a lot and it makes all calculations quite ele- 
mentary. In Section IV which is the key section of our 
Letter, we present our main lemmas and theorems for 
N— tuples of fidelities, especially, we show which form 
of the multipartite state representing QCM can be used, 
and also how one can connect fidelity calculations with 
the formalism of Young diagrams. It follows, that to de- 
termine the admissible region of fidelities, it is enough 
to consider overlaps of real vectors with the matrices 
of irreps of the symmetric group. What is important, 
all results are valid for N— tuples of fidelities, so they 
describe an asymmetric 1 — > N doner in general. In 
Section [V] which is the case study of the 1 — > 3 uni- 
versal quantum doner, we present how to obtain the 
allowed region for triples of fidelities. We also show 
that our results are consistent with: a) calculations for 
a symmetric cloning case, predicted by the Werner's 
formula 1IT21 and b) partially with results obtained in 
(161 , where optimal fidelities for this kind of LTQCM 
were presented. At the end, we give an application of 
our model, namely we show that for any given triple of 
fidelities from the allowed region, we could reconstruct 
a state that gives rise to that fidelities. This technique 
is actually quite general and could be applied to the N 
number of clones. In Section [VI] we briefly explore our 
approach with a higher number of clones, by studying 
the case of a 1 ->■ 4 UQCM. 



II. STATEMENT OF THE PROBLEM 

Let us recast a question of cloning in an equivalent 
picture of entanglement sharing (for simplicity, we fo- 
cus here on our case study-example, namely the 1 — >■ 3 
UQCM - all calculation could be easily extended to the 
case of N clones). Suppose that we have an initial state 

described by the Bell state: \xp + ) = J\ (|00) + |11)). 
We also have a cloning machine M, which could be 
described by a completely positive, trace preserving 
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(CPTP) map A (note that the cloning map A applied to 
the second subsystem of the maximally entangled state 
\ip + ), when the first is untouched, produces, in general, 
mixed state that contains all the information about the 
map). As an output, we want to obtain N shares of our 
initial state. Our scheme is presented in Figure [T] (for 
our case-study example). Let us write a state obtained 



A 



cloning 
machine M 




> Ms 



FIG. 1: Situation before and after cloning, 
after an application of the cloning map to the half of 



P1234 



(5) 



where indexes 1, 2, 3 and 4 are connected to an initial 
state and clones, according to Figure [T] 

We want to calculate an allowed region for singlet 
fractions Fu (between the initial state and one of the 
three clones (copies)), denoted by 112211 



hi 
il3 



« 2 l Tr 34(|0l234)|'/'i + 2), 

(</>13 1 Tm(puM) 

(< 4 l Tr 23(l°1234)|</>i + 4 >- 



Note that: 



hi = (ft/|Tr n (pi234)|«/'ij), 



(6) 



(7) 



where Tr^j means partial trace over all systems except 
If, and \ip^) and pn34 ar e defined below. The vector 
\xpy) = U <8> I IV 7 ^}/ is obtained after the action of 
-io-y on \xp+): 



IT) 



i(Ty\lp + 
1 



-«r y y - (|00) + |11» 



(8) 



(|01) -|10», 



where cr v is one of the Pauli matrices: cr v = 9 J 

J y [i o 

Using |8} we can write that: 

|f ~) = IT®I|?+), (9) 

where U = —iffy- The state P1234 from Eq. |7| is ob- 
tained after the following transformation: 



P1234 = (I®A)|^ lT > 



(u®i)(ot®A)|^ 1 t)^+|)(ir t ®i). 



(10) 



The four-partites states P1234, with the constraint p\ = 
1/2, are in one-to-one correspondence with cloning ma- 
chines, and the cloning fidelity of a given machine is 
determined by fidelities Fu of the corresponding state. 
Note that for any channel A, fidelity of the state 

NIg>A^ (\ip + )(ip + \) can be related to the average fi- 
delity of transmission of an initial state as follows |23J: 

f = ^ (ID 

where d is the dimension of the Hilbert space H = C rf . 
Thus instead of cloning fidelities /, we can consider 
singlet fractions F, which we will further call simply 
fidelities, while the fidelity of cloning we will term the 
'cloning fidelity'. 

Werner |12J provided the following formula for an 
optimal cloning fidelity of universal symmetric N\ — > 
N2 cloning machine: 



/n,n 2 



Ni (N 2 -Ni)(Ni + l) 
N 2 N 2 {N 1 + d) 



(12) 



In our case, when N-y = 1, N 2 = 3 and d = 2, we obtain 
that the fidelity / for the universal, symmetric cloning 
machine should be equal to 



/ 



1 4 _ 7 

3 + 9 ~ 9' 



(13) 



We can now formulate question for our case study- 
example: which values of triples of cloning fideli- 
ties (/12//13//14) are allowed for a universal cloning 
machine? (Of course in general our question is the 
following: which values of N— tuple of cloning fi- 
delities (/i2//i3/ • • -/fin) are allowed for a universal 
cloning machine?) As said above, we shall address 
equivalent question: what values of triples of fidelities 
(F12, F13, F14) are allowed for an arbitrary state of a max- 
imally mixed first subsystem? In the next sections, the 
answer is presented. 



HI. MATHEMATICAL INTRODUCTION: SCHUR - 
WEYL DECOMPOSITION 

In this section we introduce necessary mathematical 
tools from group theory. We are especially focused on 
Schur-Weyl duality ||24|. 
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Consider a unitary representation of a permutation 
group S n acting on the n— fold tensor product of com- 
plex spaces C d , so our full Hilbert space is V. = (C d )® n . 
For a fixed permutation TC £ S n a unitary transforma- 
tion V n is given by 

V n {\h} <g> . .. ® |z„)) = l^j)} (g) ... ® 1^)}, (14) 

where • is a standard basis in (C d )®". The 
space of rank— n tensors can be also consider as a rep- 
resentation space for a general linear group GL(d, C). 
Let U £ GL(d, C), thus, this induces in the tensor prod- 
uct (<C d )® n the following transformation 

U m (\h) (8) ... (8 |z n )) = LT|ti> ® .. . ® lT|z„). (15) 

A key property is that these two representations turn 
out to be each other commutants. Any operator on 
(C d ) m that commutes with all U m , Vtf £ GL(d,C), is 
a linear combination of permutation matrices V„. Con- 
versely, any operator commuting with all permutation 
matrices V n , Mn E S n , is a linear combination of U®". 
This duality is called Schur-Weyl duality. It was shown 
(see, for example, [25 1 or [26 1) that there always ex- 
ists some basis called the Schur basis which gives de- 
composition of V n and U®' 1 into irreducible representa- 
tions (irreps) simultaneously. Thanks to this, the space 
(C")® n can be decomposed into irreducible representa- 
tions of S n 

(C d )® n S 0^®7£f, (16) 

where A labels inequivalent irreps of S n and is 
the multiplicity space. It is called Schur-Weyl decom- 
position. The labels A are allowed partitions of some 
natural number n. Every partition is a sequence A = 
(Ai, . . . , A r ) satisfying 

r 

V,- A,- > 0, Ai > A 2 > . . . > A r , £A; = n, (17) 

i=i 

where r E {!,... ,n}. Every such partition corresponds 
to some diagram, which is called the Young diagram 
Il27l . Here are few examples of Young diagrams for 
n = 4, 6 and 3 respectively: 



A =(2,2), A =(3,2,1), A =(1,1,1) 

In this Letter we are interested in representations on the 
symmetric part Tif . For example the SWAP operator 
Y n can be decomposed, due to Schur - Weyl decompo- 
sition, in the following way: 

V« = 0H r(A) ®Vi (18) 

A 



where n E S„ and r(A) is the dimension of a unitary 
part. The operators are irreducible representations 
of S n . Thanks to the above-mentioned method, we can 
decompose LT®" -invariant states in the following way: 

Pl...n = 0Ir(A)®P A - (19) 
A 

Note that fidelities with singlet states are invariant un- 
der averaging over U <S> U transformations. Therefore 
the N — tuple of fidelities is invariant under an appli- 
cation of transformation U® n to the state, so we can al- 
ways use density operators that are commutant of IT®" 
i.e. they are of the form fl9) . 



IV. EXPRESSION FOR AN ALLOWED REGION OF 
N-TUPLES OF FIDELITIES. 

In this section we provide a general formula for an al- 
lowed region of N-tuples of fidelities in terms of over- 
laps of pure states with irreducible representations of 
S n (n = N + 1). This is contained in Theorem [2] We fur- 
ther show in Lemma [3] that one can restrict attention to 
pure states with real coefficients, since they determine 
the admissible region of fidelities. 

Lemma 1. Fidelity as defined in |7) is of the form 

Fik = E F i\' ( 2 °) 

A 

where 

4 = ^Tr(/ff ljt) ), (21) 

The lower index (lk) means a permutation that swaps 1 and 
k, and p^'s are arbitrary normalized states on partition A. 

Proof. From the definition of a fidelity we can write 

flic = (felpulfe) = Tr (ftfclfeXfel), (22) 



where = \{id lk - Y lk ), p lk = Tc^pi...„ and 

Tr^r denote partial trace over all systems except 1 and 
k. Expanding p2) , we obtain: 




Now we can use Schur - Weyl decomposition to repre- 
sent V( ljt) and pi.. M : 

V (1J fc) = 1,(A) ®Vf u) , Pl ... n = I r(A) ®p\ (24) 

A A 
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Inserting ( |24) into 1 23 ), we have: 
hk~- 




Equation 1 25 1 could be rewritten as: 



(25) 



(26) 



where = 2 ~ 2 Tr (i° AV (\)t))' and i° A = rf ^? A 311(1 rf A 
stands for the dimension of a given partition. 

□ 

Now we are in position to formulate the main theo- 
rem of this section: 

Theorem 2. The set T of admissible vectors of fidelities 
{F 12 , . . . , Fin} is °f the form 



T = conv 




(27) 



where conv stands for a convex hull, the union runs over all 
irreps ofS„ and 

^ A = {(f a ,...,F a ) :W6C^}, (28) 

w/iere Fj^. are of the form: F^ k = \ — 2(^^(1^^' an ^ 
where \ tp) is a pure state. 

Proof. At the beginning, let us consider the following 
mapping: 



F : P^TZ 



N 



(29) 



which maps states p\„ M 6 P(P stands for a convex set 
of all states p%,.. n ) into the N-tuples (Fi 2 , . . . , F\ n ) £ TZ N . 
Explicitly, we have 

F(pi... n ) = [F\2{p\...n), ■ ■ ■ ,F\ n {p\...n)] ■ (30) 
This mapping is affine (Lemma [9j, i.e. is of the form: 

F(pi... n ) =P(pi... n )+C / (31) 

where F : P — >■ TZ N is linear. Indeed, one can see that 
RHS of J301 could be written as: 



Fll{pi...„) = 9 Tr ( V (l)...(n) -V( 12 ))|0l...„ 



2 

1 1 



Tr (V(l2)|0l...«), 



(32) 



F\n(P\...n) 



2 Tr 
1 1 



( v (i)...(«) ~ 
Tr (V(i„)Pi...„). 



Comparing |3T 
in our case C 



and 



f32), we obtain that 
and ?(pi... n ) = 
which is ob- 



- Tr(V (12 )|0 L ..„), . . . / Tr(V( ln )pi...„) 
viously linear with respect to p\ n . 

One can note that in general |0 A from Eq. p4| | is a 
mixed state, but according to Lemma [9] we can take the 
mapping of extreme points of p\ „ ofthe form: 



= T %A) ® rf Ap A 
A flA 



0j-Ir(A)®P A = 0j 



V(A) 



(33) 



where /) A = so namely, it is normalized, and 
Tr,o A . We can now eigen-decompose p : 

p x = LpM)(ipi\, 

i=l 

where Ea Pa = 1- 

Inserting Eq. p4| into Eq. |33j, we get that: 



(34) 



0Epap a 

A z'=l 



V(A) 



(35) 



we see from | [35) that extreme points are of the form: 

J*l®|^)tyA|, (36) 



r extreme 



where A runs over all irreps from S„ and |i/? A ) is an 
arbitrary state form the space of irrep linked to a given 
partition A. Inserting Eq. |36j into Eq. |25}, we obtain 
the desired result. □ 

We note that to determine the allowed region of fi- 
delities, it is enough to consider only vectors of real 
coefficients. 

Lemma 3. To generate a convex hull ofthe allowed region of 
fidelities, it is sufficient to consider pure states of real coeffi- 
cients only. 

Proof. We need to show that, in our case, a kind of ma- 
jorization of complex pure states by real ones occurs. To 
prove that, note that our operators V A 1; j |43| are sym- 
metric and real 1 , so they could be written in a general 



1 Thanks to |28|, we know that swap operators V^ ; - ar e 
symmetric and real. It easy to see that every swap opera- 
tor between system 1 and i is of the following form: Vj x i; j = 

(Ytl2)^(23) ' ' ' ■ ' %-2,i-l)) %-lJ) {%-2,i-l) ■■■■■ ^p3)^fl2)) 

AV^. j « A T . We know that for every symmetric matrix S and an ar- 
bitrary nonzero matrix X, their composition XSX T is also symmetric. 
Thanks to this, the swap operator Vi ., is an symmetric matrix. 
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form as: 



Vlk v 2k 



Vkk 



(37) 



where k = dim'H^. Now let us write a density ma- 
trix of a pure state \ip) = («i,fl 2 , • ■ ■ / a k) T with complex 
values (letter C corresponds to the word 'complex'): 



|flx| 2 a x a 2 •■■ aia k \ 

- - i 1 9 

fl l fl 2 | fl 2P ■ ■ ■ fl 2 fl fc 



where tpj are phases connected with a given |a,-| and 
by • we denote a Hadamard product of two matrices. 
Now our goal is to show that the matrices A and C are 
positive-semidefinite. First let us define a set of vectors: 

\a)j) = cos<fy|0) +sin<fy|l), for = 1,2,..., k. (42) 

It is easy to show that matrix C can be rewritten in 
terms of vectors \u>i) i.e. Qj = (u>j\cvj) for i,j = 
1,2, ...,k. Thanks to this operation we can conclude 
that C is a square Gramian matrix. From 11291 we 
know that every square Gramian matrix is positive- 
semidefinite. 

The matrix A is of the form \(p)((p\, with |^) = 
(38) Yji KILO/ so it is positive-semidefinite. Now using 



klV 



Next, let us rewrite Tr (jj^W^^j from Lemma |lj using 
d38l and &7), as follows: 



k k 

Wj\ 2 vjj + J2 ( a j a s + a i u s) v js 

j=l j,s=l 
k k 



Tr(p£vf l0 



E \ a j\ 



7-1 



t2 E Re ( fl / fl s) 0;s- 



(39) 



We see that calculating the trace p9) with the operator 
p^. is equivalent to calculating the trace 1 39 1 with the 
following matrix p^-c 



Pr-c 



\a\\ 2 Re(fli«2) ■■■ Re(«iflfc) 
Re(fli«2) | fl 2| 2 ••■ Re( fl 2 fl fc) 



Re(fliflfc) Re(fl2«fc) ■■■ |«jt| 



(40) 



where the new index R — C means that we are left only 
with real numbers. We shall now show, that p^_ c is 
positive-semidefinite. 

Let us denote the total phase in this case by e^i for 
j — 1, . . . , k, so we are left with Re(e 1 ^' ) = Re(cos <pj + 
i sin (pj) = cos <pj. We then obtain 



Pr-C = 



|fl 1 fl 2 | \a 2 \ 2 



|fl 2 «fc| 



fli«)tl Mjtl •■■ kl 2 

1 C0S((^1 - <£ 2 ) • • ■ COS(01 - ^Jt) 

COs(<^!-^2) 1 ... COs((p 2 -<pk) 



C0S(^! -tf>, f ) COs((p 2 ~(pk) 
~-A*C, 



1 



(41) 



Fact 
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we obtain that our operator p'\_ c is also 
positive-semidefinite. 

Because p y \_ c is real and symmetric we get from 
Fact [ll] that the matrix p^-c possesses real eigenvec- 



tors, so indeed it is a mixture of real pure states. 



□ 



V. CASE STUDY: A REGION FOR TRIPLES OF 
FIDELITIES 

A. Partitions and transpositions 

In our case-study example we have four particles 
(n = 4) which means that allowed partitions are Ai = 
(4),A 2 = (3,1),A 3 = (2,2),A 4 = (2,1,1) and A 5 = 
(1,1,1,1). Because we want to consider only qubits 
here, so in our case d = 2, and hence A runs over binary 
partitions only or, equivalently, over Young diagrams 
with two rows (for other diagrams the multiplicity 
space becomes zero-dimensional). Thanks to this, 
we are left only with A a = (4),A 2 = (3, 1),A 3 = (2,2). 
For the partition A2 = (3,1) unitary representations of 
transpositions T(l,2), T(l,3), T(l,4), which we shall 
need, are |30|: 



j\ 2 

(12) 



rA 2 

(14) 



1 
1 
0-1 



^2 

(13) 



1 

U 2 2 

1 

U 2 2 . 



_1 _s/2 _\/6' 

3 3 3 

V2 5 _ \/3 

3 6 6 

\/6 _V3 1 

3 6 2 ■ 



(43) 



For the partition A3 = (2, 2) unitary transposition repre- 
sentations T(l, 2), T(l,3), T(l,4), which we shall need, 
are 13!: 



'A3 

(12) 
(14) 



1 
-1 



_1 ^3 

2 2 

^3 1 

2 2 J 



f^3 
(13) 



_ 1 _\/3 
2 2 
\/3 1 
2 2 . 



(44) 
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Note that for the partition Ai we have a trivial represen- 
tation, so it is not reported here, explanation is provided 
later. 



B. Partial result: fidelity regions for each partition 



According to Lemma [9] all extreme points are just 
the convex hull of the figure from Figure [2] It al- 
lows us to use pure states only. In our situation, we 

use states which are of the form |i/>( 2 ' 2 ') - 1 



«2 



and 



£ TZ, so they are real pure states (see 



Lemma |3|. Each of this state generates a pure state p 



of the form p 



(2,2) 



1 

a\U2 



, where a 2 



1, for 



the partition A = (2, 2), and p 



(3,1) 



0-10-2 
«1«3 



0\02 0\Oo, 

„2 



(7, 



where a\-\- a\ + a\ = 1, for A = (3, 1) respectively. In- 
serting p( 2 ' 2 ) and pt 3 ' 1 ) into Equation pO} from Lemma 
[l] we obtain the plot 2 in Figure |2j with values F\2, F13 
and F14 on axes. Now we want to give an explanation, 
why the partition A = (4) is not included in figures. It 
can be shown that in this case, all fidelities are equal to 
0, and since we do not want to 'blur' the main message 
of this Letter, this partition is excluded from Figure \2 
But, strictly speaking, in the next step, the convex hul 
also with this point should be created. This should lead 
to a 'larger' possible region for fidelities. Of course, it 
is not a hard task from the computational point of view, 
since taking a convex hull of some point and some fig- 
ure is quite an elementary operation. 
Remark 4 Note that this does not mean that cloning fideli- 
ties are equal to and then that the partition A = (4) should 
he excluded because perfect universal anti-cloning is forbid- 
den by the laws of quantum mechanics. In our case singlet 
fractions are equal to and fidelities are equal g (see Eq. ( |TT) 
and the comment in the paragraph below it). 

At the end of this paragraph we present also explicit 
formulas for fidelities for every irrep labeled by parti- 
tions A2 and A3. 





FIG. 2: Plot of Fu = Ea F ii with i° (2,2) and (° (3,1) • On axes val- 
ues of F12, F13 and F14 are reported. The ellipse corresponds to 
the partition A = (2, 2), when the other figure to the partition 
(3,1). On the ellipse, a small numerical error can be seen. 
On each plot three coordinates are presented from the set: 
A = (0,1,0), B = (1,1,0) and C = (0,0,0), and C' = (1,1,1). 
Two different views of plot are presented. 



2 A11 plots are obtained using Mathematica software. 
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C. The main result: an allowed region for fidelities 

Since, we have two partitions (2, 2) and (3, 1) and cor- 
responding pure states p^f^ 1 and , we obtain two 
figures in Figure [2] But, we are interested in situation, 
where we can obtain a general answer to our question 
from Section |TTJ namely in situation where we have a 
mixture of both partitions: £a Pa^iV To solve this task, 
we can construct a convex hull of figures from Figure [2] 
It is presented in Figure [3] One can also see that our Fig- 
ures |2] and [3] are invariant under rotation of angle 2n/3 
along straight line Fj2 = £13 = Fi4- This corresponds to 
three conjugacy classes 3 . 

Remark 5 Because of the properties of the cloning map A 
(see Sec. all possible convex mixtures of the partitions 
shown in Fig. ^are possible and produce a correct quantum 
doner, i.e. a trace preserving completely positive map. 




D. Comparison with other models 

As we mentioned before, this particular UQCM - an 
asymmetric 1 — >■ 3 doner - was studied before in lTT6l 
and in (17l . It was shown there how to obtain an opti- 
mal region for this particular doner. Plotting Eq. (38) 
for d = 2 together with a constraint from Eq. (32) from 
[16 1 and then using Eq. (H\ (from this work) to obtain 
fidelities, one can check that then the optimal doner 
corresponds to a region for fidelities obtained for the 




FIG. 3: The convex hull of figures that corresponds to the 
partition (2, 2) and (3, 1), based on Figure[2] On axes values of 
F12, F13 and F14 are reported. On each plot, three coordinates 
are presented from the set: A = (0, 1, 0), B = (1, 1, 0) and C = 
(0,0,0), and C = (1,1,1). Two different views are presented. 



3 Two permutations re and n are conjugate iff n = C • n ■ <T 1 , where 
u is also a permutation and a, n, n 6 S„ . 
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partition A = (3, 1 ) (but first some unknown at this mo- 
ment constraint should be applied, leading to a cut-off 
of the plot of the partition A = (3,1), so it will corre- 
spond to the result from |16|). It is an important remark 
that one should keep in mind - in our work optimal 
doners are not considered, instead, the whole possible 
region for a given N— tuple of fidelities is obtained that 
later somehow could be optimized by 'throwing away' 
some partitions and establishing some constrains. What 
is more, authors of the paper [16J also proved the opti- 
mally of the asymmetric 1 — > 2 quantum doner using 
their technique. This result is also in accordance with 
our results, one can check that using our method, by 
taking valid partitions from the symmetric group S(3), 
the same result could be obtained. In general, expres- 
sions for possible fidelities in both cases (so our and 
from J 16J) look quite similar (compare for example 'our' 
Eq. (g6) with Eq. (38) from ESI), but because of the 
constraint from Eq. (32), the optimal fidelities could 
be obtained, not the full region like in our work, so it 
would be worthy to explore the direct connection be- 
tween these two approaches in the future. 



E. Special case: symmetric cloning 

It can also be shown that in the case of symmet- 
ric cloning, the maximal possible value of the triple 
{Fi2, £13, F14) is (5/5/3)/ so it is in accordance with 
the formula obtained by Keyl et al. |13| and Werner 
llT2l (first using Equation dTTJ), namely we are able 
to find such a point (F12, F13, F14) that corresponds to 
the case of the optimal symmetric fidelity F. To find 
the maximal possible fidelity (for a triple (F 12 , F i3 , F 14 )) 
which corresponds to the case of the optimal symmetric 
cloning, one method is to find a plane equation which 
includes the ellipse (which corresponds to the partition 
(2,2)). It can be shown that it is given by the formula 
F12 + F13 + Fu = §- It can be obtained, for example, by 
taking coordinates of points from Figure [2] and then cal- 
culating the plane equation according to the following 
well-known formula: 



x-xi y-yi z-z\ 

X 2 -Xi 1/2 - J/l Z 2 - Zi 
*3 -X 2 J/3 - J/1 Z 3 - Zi 



0. 



(47) 



So, of course we have to identify x, y, z with 
^12/ ^13 an d F14 first, and then insert to | [47) val- 
ues of three points Fi(F 12 ,F 13 ,F 14 ), P 2 (F 12 ,F i 3 ,F 14 ) and 
^M-FU/ F13, F14) from the ellipse and solve equation 
which is then obtained. 



F. Applications 

The above-mentioned approach allows us to recon- 
struct states from every subspace Hf which satisfies 



not only the symmetric case (Fj = F2 = F3) (like in 
a typical QCM) but also more general relations be- 
tween fidelities. In this section we reconstruct states 
belonging to the subspaces corresponding to partitions 
Ai = (2, 2), A2 = (3,1) and satisfying the following for- 
mula 



max (Fj + F 3 = 2F 2 ) . 



(48) 



a) Partition A x = (2,2) 

In this case we have two constrains Fj + F3 = 2F2 and 
a 2 + «2 = 1, so after solving this system of equations we 
obtain four allowed pairs of solutions 

(«1,«2)(1) = 



(49) 




Fidelities corresponding to the above pairs are 



( r-A-l Ml F^i \ 
l^l ,t 2 ,t 3 J (1) 

( pAi pAj F Aj \ 
\h > b 2 ' b 3 )(2) 

(^2^3%) 

(FM Al ^3 Al )(4) 



1(2-V3),i,i(2 + V3) 

1(2 + V3),i,|(2-V3) 
\ (2 -V3), \,\(2 W3) 



^(2 + V3),i,i(2-V3)). 



(50) 



We can see that the maximal fidelity F 1 1 can be obtained 
for a pair (2) and (4), so our reconstructed states are the 
following 



Ai 
P(2) 



1/2-V3 -1 
4^ -1 2 + V3J' 

1(2-V3 1 \ 
4V 1 2 + V3J' 



(51) 



Now note something more general. Thanks to Lemma 
[lj we have F^ 1 = a\ and then a 2 = 1 — a\ = 1 — F^ 1 , so 
we can express every states jO Al 6 Wf in terms of the 

fidelity Fj Al 



1 - F Al 



Ff 1 (1 



f Ai ; 



±VF Ai a-fr) 



(52) 
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For example, the sign '+' in p2) corresponds to states 
obtained from pairs (a±,b£) and (^3,^3) while the sign 
'— ' corresponds to states obtain from pairs (a\,b\), 
(«2/ h)- 

b) Partition A 2 = (3,1) 

For the partition A 2 = (3,1) we have a more compli- 
cated situation. Here we have three parameters a\, ai, a$ 
but only two constrains F1+F3 = 2F 2 and «? + «| + °3 
1 , so we first need to eliminate the parameter a\ and 
then express the parameter a 3 as a function of a 2 . Nu- 
merically, we find that the fidelity F\ is maximal for the 
following values of parameters fli,«2/ fl 3 



need, are |30 



(«l/«2/«3)(i) = (0.114,0.318,0.941), 
(«i,«2,«3)( 2 ) = (-0.114,-0.318,-0.941), 
(fli,fl 2 ,fl 3 )(3) = (0.114,-0.318,-0.941), 
(fll,«2,«3)(4) = (-0.114,0.318,0.941). 

Corresponding states are 



(53) 



V (12) 



V (14) 



yA 2 
V (15) 



"1 








- 





1 














1 











- 


-1 


"1 














1 

3 









V ^ 

3 









vl 








1 


1 






4 


4 




1 


/5 






4 


V 3 








y§ 






2 








1 








2 


V 2 





fX 2 
(13) 



x/2 
3 

5 
6 

1 

2V3 

5 
3 

U 
12 

1 

6^2 
1_ 
2V6 



-1 








- 





1 














1 

2 


V3 
2 


.0 





73 
2 


1 

2 J 



_ 

2 
3 
1 

2x/3 

1 

2 

1 /s 

2 V 6 
1 

6\/2 

5 

6 

1 

2a/3 



1 5 

2 V 2 

1 

2v/6 
1_ 
2v/3 
1 
2 



(55) 



A2 A2 



A 2 A 2 
P(3) P(4) 



'0.013 0.036 0.107N 
0.036 0.101 0.299 
,0.107 0.299 0.886, 

0.013 -0.036 -0.107N 
-0.036 0.101 0.299 
-0.107 0.299 0.886 



(54) 



Note that fidelities in this case are IF 



^2 pA 2 pA 2 



-! ' x 3 

(0.886,0.556,0.220). The numbers for our numerical 
data are obtained by Wolfram Mathematica Alpha |31 1. 
Most likely they are indeed optimal. 



Now, assuming that our pure state is of the form 
l^/ 4 - 1 )) = [fl 1 ,a 2 ,fl3,fl4] T , we obtain from Lemma [l] the 
following formulas for fidelities: 



12 



13 



P \ 2 
14 



K 1 



3 ' 4 / ' 

^ + V3fl 3 fl 4 - y 



VI. FIDELITIES FORMULAS FOR A 1 
QUANTUM CLONER 



4 UNIVERSAL 



In this section we briefly report formulas for fideli- 
ties of a 1 — »■ 4 LfQCM for every allowed irrep A for 
an S(5) symmetric group in the case of qubits. In 
this particular example we are left with three parti- 
tions only i.e. Aj = (5), A 2 = (4,1) and A 3 = (3,2). 
For the partition A 2 unitary representations of trans- 
positions T(1,2),T(1,3),T(1,4),T(1,5), which we shall 



pA 2 
r 15 



«3^4 

" V3 



«2«3 
3^2 



.-«1«3+ 




fl 2 fl4 ^3^4 



(56) 



with a normalization condition fl? = 1. 
For the partition A3 unitary representations of trans- 
positions T(1,2),T(1,3),T(1,4),T(1,5), which we shall 
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need, are |30|: 



7^3 
(12) 



V 



(13) 



"1 





o - 






n 
U 


1 U U 


n 
U 









10 












1 















-1_ 






-1 











o - 







_1 

2 


_V3 
2 













V3 
2 


1 

2 



















1 

2 


V3 
2 




.0 








73 

2 


1 

2 - 






1 

3 


3 











V 3 = 

V (14) 



x/2 
3 



5 
6 

1 

2\/3 



(15) 



















1 


1 


1 


3 

1 


3 a/2 
1 


^6 
1 


3\/2 


6 

1 


f 


76 

1 


v'3 
1 


2 
1 


x/6 
1 

sfl 


75 

1 

2 


2 





1 

2\/3 

1 

2 






_ j_ 

"V6 
_ j_ 

~f 

2 
1 
2 









_i a/3 

2 2 
\/3 1 
2 2 . 

j_' 

V2 



(57) 



Assuming that our pure state in this partition has 

a form li/>( 3 ' 2 )} = [ai,fl2/ fl 3/ fl 4/ fl 5] T we obtain from 
Lemma u\ the following formulas for fidelities: 




5 ' 



+ V3, 



-V3< 



«4«5 - y 



«2«3 



2\/2flifl2 

6~ 




-2W ^flifl3+ 



2«2«3 



3 ^ 4 + ^_ 



V3 



+a 2 a 5 



(58) 



with a normalization condition Ef=i fl f — !• After tak- 
ing the convex hull, the allowed region for fidelities 



could be obtained also for this quantum doner. Of 
course, the state reconstruction technique from Section 
IV F| could also be used in this case. At the end of this 
section, let us note that here we also do not consider 
partition Aj (see Section VB| . 



VII. CONCLUSIONS 

We have shown that by using representation theory, 
especially Young diagrams, action of the universal 1 — > 
N quantum cloning machine could be described. The 
method of irreps is quite powerful, because it allows 
us to decompose (usually big) Hilbert space into blocks 
(linked with a given partition A) of smaller dimensions 
which, of course, are easier to deal with. For example, 
in our case-study example 1 — > 3 UQCM, the Hilbert 
space C 16 is decomposed into blocks of dimension 1, 
2 and 3 respectively. What is more, using our model, 
fidelity expressions are quite easy to obtain, one only 
needs to know representations of all possible irreps for 
a given symmetric group. 

We have also shown that the convex hull could be 
made to obtain full knowledge about our model, in 
the sense that it gives raise to the full possible range 
of fidelities. After careful studies, also optimal UQCM 
could be found. What is more, we have shown that 
by restricting ones attention only to real pure states in 
each of the block, the full answer to our initial question 
is obtained. We point out that our UQCM gives correct 
results for the case of the symmetric cloning and we 
have also proved that the optimal value, in the symmet- 
ric case, could be obtained. What is more important, 
our approach allows to reconstruct any given state con- 
nected to a given N— tuple of fidelities. 

Of course, in future, it would be interesting to ex- 
tend our method in such a way that also qudits could 
be described. What is more, it would be interesting to 
search for method that allows to obtain the optimal N- 
tuple of fidelities - comparison of our results with these 
obtained in [16J may suggest that the optimal fidelities 
always correspond to only one partition (with some cut- 
offs) A, since for cases 1 — > 2 and 1 — > 3 we have that 
in the first case, the optimal fidelity region corresponds 
to A = (2,1), in the second one to A = (3,1). Also 
it is worthy to explore the direct relation between these 
techniques, so namely 'our' technique and this reported 
in EH. 

Note added: After the completion of this paper we be- 
came aware that similar results for a 1 — » 3 universal 
quantum doner have been reported in Il32l . 
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VIII. APPENDIX 

To prove our results we need the following defini- 
tions and lemmas: 



Definition 6. A set Z is convex iff 
V X/ y g z V <,,<i yX + (1 - ji)Y e z. 

Definition 7. L(-) : x — > y is affine map iff 3 a M x L{x) 
L(x) + a, where L is linear. 

Now let us present two useful lemmas 
Lemma 8. Suppose that y, y 6 L(fi), then: 

ay + (1 — a)y £ L(fi) 
(it means that L(Q) is a convex set). 
Proof. Let us write: 

L ( x ) = y = l (Ep' x ') = YLv&ixi) + «/ 

i i 



(59) 



Fact 10. T/ze Hadamard product of two positive-definite ma- 
trices is again positive-definite f&M . 

Fact 11. If A is a real symmetric matrix, then all of its eigen- 
values are real, and eigenvectors can always be chosen to be 
real. 



Proof. Here, we only prove the second part of Fact 11 
since the first one is obvious. Let A be a real eigenvalue 
of the matrix A with an associated eigenvector v. Let 
us write v as v = Re(i>) + flm(c), so that A(Re(v) + 
Im(iO) = A(Re(p) + Jm(v)). It implies that ARe(v) = 
ARe(v) and Alm(v) = \lm(v). Now, if Re(v) ^ 0, 
then it is a real eigenvector of A. When Re(c) = 0, then 
Im(c) is a real eigenvector. □ 



(60) 



where we use the fact that L is affine and x = Y^i Pi%i, 
x — JZiqiX for some x i ,x' i 6 E(O). Inserting |60j into 
( |59) we get: 

oi(Y^PiL(xi) + a) + (1-«)(X^/L (*,•) + «) 

i i 

= aa + (1 — a)a + a^L(p,x,) + (1 — a) ^L(^,x ; ) 

i i 

= a + L(ocx + (1 - a)x) = L(ax + (1 - a)x') G L(Q). 

(61) 



This ends the proof of lemma. 



□ 



Lemma 9. Suppose that we have an affine map L : X — > Y. 
If by Ci we denote a convex subset of X, where X is a finite 
dimensional space and by E(O) a set of extreme points ofCl, 
then L(n) is a convex set and L(E(tl)) can reproduce set 
L(Q) after taking the convex hull, i.e.: 



L(n) = conv L(E(Q)). 



(62) 



Proof. We can write L(E(Q)) C L(Q), because we know 
that E(fi) C n. This together with Lemma ^ ends our 
proof. □ 
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